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A  COMPUTATIONAL  STUDY  OF  THE  STRUCTURES  AND 
ELECTROSTATIC  POTENTIALS  OF  SOME  AZINES  AND 
NITROAZINES 


JANE  S.  MURRAY,  JORGE  M.  SEMINARIO  and  PETER  POLITZER 
Department  of  Chemistry,  University  of  New  Orleans,  New  Orleans,  LA  70148  (U.S.A.) 
(Received  22  April  1988) 


ABSTRACT 

An  ab  initio  SCF  computational  approach  is  used  to  study  six  azines  ( azabenzenes )  and  their 
mononitro  derivatives.  Our  primary  interest  is  in  determining  how  the  reactive  properties  of  the 
azines  are  affected  by  the  introduction  of  the  nitro  group.  All  structures  are  optimized  at  the  3- 
21G  level,  and  these  are  then  used  to  compute  the  STO-5G  molecular  electrostatic  potentials. 
Among  the  various  isomers,  those  having  two  adjacent  ring  nitrogens  are  by  far  the  least  stable. 
The  nitro  derivatives  are  most  stable  when  the  substituent  is  (3  to  a  ring  nitrogen  and  least  when 
it  is  a.  The  dominant  features  of  the  electrostatic  potentials  of  these  molecules  are  the  large  and 
strong  negative  regions,  centered  in  the  molecular  planes,  that  are  associated  with  the  ring  nitro¬ 
gens  and  are  indicative  of  their  basic  characters.  These  negative  potentials,  and  correspondingly 
the  basicities,  become  weaker  as  the  number  of  ring  nitrogens  increases  and  also  with  the  substi¬ 
tution  of  the  electron-withdrawing  nitro  group.  The  regions  above  and  below  the  rings  become 
increasingly  positive  in  going  from  the  mono-  to  the  tetra-azine  and  with  the  introduction  of  the 
-NO-,  group,  suggesting  enhanced  susceptibility  to  nucleophilic  attack. 


INTRODUCTION 

The  azines,  or  azabenzenes,  are  an  interesting  class  of  molecules  which  can 
be  regarded  as  related  to  benzene  through  the  replacement  of  one  or  more 
carbons  by  nitrogens.  The  presence  of  the  latter  introduces  a  degree  of  basicity, 
but  the  electron -attracting  tendencies  of  the  nitrogens  reduce  the  aromaticity 
of  the  ring  and  its  reactivity  toward  electrophiles;  indeed  there  develops  a  sus¬ 
ceptibility  toward  nucleophilic  attack  [  1  ] . 

Our  primary  interest  in  this  work  is  in  how  these  properties  are  affected  by 
the  introduction  of  the  nitro  group,  which  is  very  strongly  electron-withdraw¬ 
ing,  primarily  through  induction  [2].  Accordingly  a  computational  analysis  of 
the  azines  I-VI  and  their  mononitro  derivatives  VII-XVII  has  been  carried 
out.  For  comparison,  nitrobenzene,  XVIII  was  also  included.  For  all  these  mol¬ 
ecules,  structures  and  electrostatic  potentials  have  been  computed,  the  latter 
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being  well-established  as  a  means  for  interpreting  and  predicting  molecular 
reactive  behaviour  [3-5]. 


I  ,  pyridine 


II  ,  pyridazine 


III,  pyrazine  IE,  pyrimidine 


E  ,  s-criazine  EE,  s-cetrazine 


EH 


METHODS 


An  ab  initio  self-consistent-field  molecular  orbital  procedure  (GAUSSIAN 
82  [6] )  was  used  to  determine  the  optimized  structures  of  molecules  I-XVIII. 
The  optimizations  were  carried  out  at  the  3-21G  level,  which  has  been  shown 
to  be  effective  for  this  purpose  [7] . 

The  electrostatic  potential  V(r)  that  is  created  in  the  space  around  a  mole¬ 
cule  by  its  nuclei  and  electrons  is  expressed  rigorously  by  eqn.  (1 ) 


V(r)  =  S 


a  I Ra -r I 


0(r'  )dr' 

J  lr'-r| 


(1) 


where  ZA  is  the  charge  on  nucleus  A,  located  at  i?A ,  and  p(r)  is  the  electronic 
density  function,  computed  from  the  molecular  wave  function. 

The  electrostatic  potential  has  been  demonstrated  to  be  a  useful  tool  for 
gaining  insight  into  the  reactive  properties  of  molecules,  especially  toward 
electrophilic  attack  [3-5].  The  first  term  on  the  right  side  of  eqn.  (1)  is  the 
positive  contribution  of  the  nuclei,  while  the  second  represents  the  negative 
effect  of  the  electrons.  Thus  an  approaching  electrophile  will  tend  initially  to 
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go  to  those  regions  in  which  V7(r)  is  negative,  and  particularly  to  the  points  at 
which  its  most  negative  values  (local  minima)  occur.  It  is  important  to  note 
that  the  electrostatic  potential  is  a  real  physical  property,  which  can  be  deter¬ 
mined  experimentally  as  well  as  computationally  [5]. 

V{r)  was  calculated  for  molecules  I-XVIII  in  terms  of  the  STO-5G  basis 
set,  using  the  3-21G  optimized  geometries.  The  former  basis  was  found  to  be 
more  reliable  for  properties  related  to  the  electronic  charge  distribution. 


N  _  N 


Distances  (A) 


a:  1.331  (1.338) 
b:  1.383  (1.394) 
c:  1.384  (1.392) 


a:  1.355  (1.330) 
b:  1.316  (1.341) 
c:  1.395  (1.393) 
d:  1.365  (1.375) 


a:  1.331  (1.339) 
b:  1.381  (1.403) 


a:  1.329  (1.340) 
c:  1.332  (1.340) 
d:  1.382  (1.393) 


a:  1.330  (1.319) 


a:  1.329  (1.334) 
b:  1.332  (1.321) 


Angles  (deg) 


e-a:  119  (116.9) 
a-b:  123  (123.8) 
b-c:  119  (118.5) 
c-d:  119  (118.4) 

f-a:  119  (119.3) 
b-c:  123  (123.7) 
c-d:  117  (117.1) 


c-a:  118  (115.6) 
a-b:  121  (122.2) 


a-b:  125  (127.6) 
b-c:  118  (115.5) 
c-d:  122  (122.3) 
d-e:  117  (116.8) 

a-b:  124  (126.8) 
b-c:  116  (113.2) 


c-a:  125  (127.4) 
a-b:  118  (116.0) 


F.g.  i.  Calculated  strictures  of  the  unsubstituted  azines  I- VI.  Experimentally  determined  values 
are  in  parentheses  ( 1  ].  All  the  molecules  are  planar. 
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3zn 


Discances  (A) 


a : 

1.305 

f : 

1.329 

b: 

1.374 

g: 

1.460 

c : 

1.383 

h: 

1.228 

d: 

1.384 

i: 

1.253 

e : 

1.385 

Angles 

(deg) 

a-b:  124 

f-a : 

119 

b-c:  118 

b-g: 

119 

c-d:  119 

h-i: 

126 

d-e:  118 

g-h: 

118 

e-f:  119 

xrn 


IX 


X 


XE 


a : 

1.330 

h-a : 

119 

b: 

1.382 

a-b : 

122 

c : 

1.375 

b-c : 

117 

d: 

1.454 

c-g: 

121 

e : 

1.242 

c-d : 

119 

e-f  : 

126 

a: 

1.326 

f : 

1.333 

a-b : 

121 

f-a: 

119 

b: 

1.379 

g: 

1.440 

b-c : 

120 

b-g: 

120 

c : 

1.378 

h: 

1.243 

c-d : 

118 

g-h: 

117 

d: 

1.380 

i: 

1.244 

d-e : 

118 

h-i: 

125 

e : 

1.383 

e-f : 

123 

a: 

1.350 

f:  1.319 

a-b : 

119 

f-a:  119 

b: 

1.294 

g:  1.462 

b-c: 

125 

b-g':  117 

c : 

1.387 

h:  1.226 

c-d: 

116 

g-h:  118 

d: 

1.364 

i:  1.253 

d-e : 

118 

h-i:  126 

e : 

1.395 

e-f : 

123 

a: 

1.356 

f : 

1.315 

a-b : 

120 

f-a:  120 

b: 

1.314 

g: 

1.444 

b-c : 

122 

c-g:  120 

c : 

1.388 

h: 

1.241 

c-d : 

119 

g-h:  117 

d: 

1.357 

i : 

1.242 

d-e : 

116 

h-i:  126 

e : 

1.395 

e-f : 

123 

Fig.  2.  Calculated  structures  of  the  nitroazines  VII-XI.  All  the  molecules  are  planar. 


RESULTS  AND  DISCUSSION 


Structures  and  stabilities 

The  key  features  of  the  3-2 1G  optimized  structures  of  the  azines  and  their 
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Distances  (A)  Angles  (deg) 


.306 

f  : 

1.330 

a-b : 

123 

t-a : 

118 

..377 

g: 

1.452 

b-c : 

119 

a-g : 

117 

..329 

h: 

1.227 

c-d : 

119 

g-h: 

118 

..332 

i: 

1.252 

d-e : 

121 

h-i : 

126 

L.382 

e-f  : 

120 

a : 

1.306 

h-a : 

127 

b: 

1.332 

a-b : 

117 

c : 

1.376 

b-c : 

121 

d: 

1.470 

c-g: 

118 

e : 

1.240 

h-d: 

117 

e-f : 

126 

328 

f:  1.333 

a-b : 

124 

f-a : 

118 

330 

g:  1.464 

b-c : 

118 

c-g: 

117 

303 

h:  1.251 

c-d: 

123 

g-h: 

116 

372 

i:  1.226 

d-e : 

115 

h-i: 

126 

383 

e-f : 

122 

a:  1.331 

a-b : 

125 

b:  1.331 

b-c : 

118 

c:  1.327 

c-d: 

120 

d:  1.378 

d-e : 

119 

£:  1.432 

d-f : 

121 

g:  1.243 

g-h: 

126 

Fig.  3.  Calculated  structures  of  the  nitroazines  XII-XV.  All  the  molecules  are  planar. 


nitro  derivatives  are  shown  in  Figs.  1-4.  Their  total  energies  and  the  relative 
energies  within  the  various  groups  of  isomers  are  listed  in  Table  1. 

For  the  unsubstituted  azines  the  experimentally  determined  structural  data 
are  available  [1],  and  are  also  given  in  Fig.  1.  Our  calculated  results  are  in 
generally  good  agreement  with  these,  the  largest  differences  being  0.025  A  and 
2.8°.  The  C-C,  C-N  and  N-N  bond  lengths  are  intermediate  between  their 
typical  single-  and  double-bond  values  [8],  consistent  with  these  molecules 
having  some  degrees  of  aromatic  character  [1  ];  in  particular,  most  of  the  C-C 
distances  are  very  similar  to  that  in  benzene.  With  regard  to  the  nitro  deriva¬ 
tives,  the  principal  structural  effect  of  introducing  the  -N02  group  is  to  shorten 
the  bonds  leading  to  the  site  of  substitution. 
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Distances  (A) 


Angles  (deg) 


a : 

1.308 

f  : 

1.333 

a-b:  126 

f-a : 

116 

b: 

1.307 

g: 

1.469 

b-c :  116 

a-g: 

117 

c : 

1.334 

h: 

1.235 

c-d:  123 

h-i: 

127 

d: 

1.328 

d-e:  117 

e : 

1.329 

e-f:  123 

a: 

1.312 

f-a : 

126 

b: 

1.329 

a-b : 

117 

c : 

1.312 

b-c : 

117 

g: 

1.459 

c-d : 

125 

h: 

1.235 

f-g: 

117 

h-i: 

127 

a : 

1.378 

e-a : 

122 

b: 

1.381 

a-b : 

119 

c: 

1..386 

b-c : 

120 

f : 

1.449 

c-d : 

120 

g: 

1.244 

e-f : 

119 

g-h: 

125 

Fig.  4.  Calculated  structures  of  the  nitroazines  XVI  and  XVII,  and  of  nitrobenzene,  XVIII.  All 
the  molecules  are  planar. 


The  relative  stabilities  of  the  various  isomers  show  certain  definite  patterns 
(Table  1).  It  is  clear  that  having  two  nitrogens  adjacent  to  each  other  in  the 
ring  confers  considerable  instability  upon  the  system;  II,  X  and  XI  are  all  much 
less  stable  than  their  isomers.  For  the  unsubstituted  diazines  II-IV,  this  has 
earlier  been  observed  both  experimentally  [9]  and  computationally  [10,11], 
and  the  same  conclusions  were  reached  in  a  recent  study  of  the  1,4-dioxadi- 
azines  and  their  mononitro  derivatives  [12].  As  before,  we  suggest  that  an 
important  factor  in  this  instability  may  be  a  tendency  of  these  systems  to  de¬ 
compose  through  the  evolution  of  N2. 

Substitution  of  a  nitro  group  does  not  change  the  order  of  stabilities  of  the 
diazines,  which  remains  II  <  III  <  IV.  However,  certain  preferences  for  the  site 
of  substitution  can  be  discerned  for  both  the  nitropyridines  and  the  nitrodi- 
azines.  The  derivatives  are  most  stable  when  the  nitro  group  is  on  a  carbon  (3 
to  a  ring  nitrogen,  and  least  when  on  the  a  carbon.  This  may  reflect  the  elec¬ 
tron-attracting  tendency  of  the  nitrogen,  which  depletes  the  electronic  density 
around  the  a-  and  y-positions  (see  structures  XIX-XXI,  for  example)  and 
makes  them  less  favorable  sites  for  the  electrophilic  nitro  group.  However,  the 
effects  of  these  site  preferences  upon  the  stabilities  of  the  nitrodiazines  are 
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TABLE  1 


Calculated  total  and  relative  energies3 


Molecule 

Total  energy,  3-21G 
(hartree) 

Relative  energies  of  isomers 
(kcal  mol-1)8 

I 

-245.3120 

II 

-261.1597 

29.2 

III 

-261.1975 

5.5 

IV 

-261.2062 

0.0 

V 

-277.1011 

VI 

-292.8894 

VII 

-447.5980 

1.8 

VIII 

-447.5982 

1.7 

IX 

-447.6009 

0.0 

X 

-463.4394 

32.3 

XI 

-463.4414 

31.1 

XII 

-463.4792 

7.3 

XIII 

-463.4811 

6.1 

XIV 

-463.4865 

2.8 

XV 

-463.4909 

0.0 

XVI 

-479.3703 

XVII 

-495.1498 

"Within  each  of  the  three  groups  of  isomers  (II-IV,  VII-IX,  and  X-XV),  the  relative  energies 
are  given  in  terms  of  the  most  stable  isomer  being  assigned  a  value  of  zero. 


much  weaker  than  the  striking  destabilization  due  to  two  adjacent  ring  nitro¬ 
gens  (compare  the  energies  of  X  and  XI  to  those  of  XII-XV  in  Table  1). 


•»5, 


IXX.  3S  TXX 

Electrostatic  potentials  of  unsubstituted  azines 


+ 
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The  electrostatic  potentials  of  the  unsubstituted  azines,  I- VI,  have  been 
presented  and  discussed  previously  [  10 }  and,  therefore,  only  certain  key  points 
will  be  mentioned  here.  The  dominant  features  are  the  large  and  strong  nega¬ 
tive  regions  centered  in  the  molecular  plane,  that  are  associated  with  the  ring 
nitrogens  and  can  be  attributed  to  their  lone  pairs.  (See  Fig.  5,  for  example.) 
Due  to  the  electron-withdrawing  power  of  these  nitrogens,  nothing  remains  of 
the  negative  potentials  that  are  found  above  and  below  the  ring  in  benzene 
[  13  ]  ( except  to  a  small  extent  in  the  case  of  pyridine,  I ) .  These  regions  become 
increasingly  positive  in  going  from  I  to  VI.  This  is  consistent  with  observations 


■  -84.1 


■  -70.2 


•  -30.1 


♦  -48.0 


Fig.  5.  Electrostatic  potential  (in  kcal  mol-1 )  in  the  molecular  plane  of  pyrimidine,  IV.  Dashed 
contours  correspond  to  negative  values.  The  positions  of  the  most  negative  potentials  are  indi¬ 
cated,  and  the  magnitudes  are  given  below  the  figure. 

Fig.  6.  Electrostic  potential  (in  kcal  mol-1)  in  the  molecular  plane  of  5-nitropyrimidine,  XV. 
Dashed  contours  correspond  to  negative  values.  The  positions  of  the  most  negative  potentials  are 
indicated,  and  the  magnitudes  are  given  below  the  figure. 


if  -50.1  O  •S0'6  0  '30-1 

Fig.  7.  Electrostatic  potential  (in  kcal  mol-1)  in  the  molecular  plane  of  3-nitropyridazine,  X. 
Dashed  contours  correspond  to  negative  values.  The  positions  of  the  most  negative  potentials  are 
indicated,  and  the  magnitudes  are  given  below  the  figure. 
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TABLE  2 


Some  properties  of  unsubstituted  azines8 


Molecule 

Electrostatic  potential  minimum 
(kcal  mol-1 ) 

P  Ka 

© 

I 

-92.9 

5.2 

© 

n 

-90.1 

2.1 

m 

-84.1 

1.1 

n z 

-82.4 

0.4 

o 

N<>N 

3T 

-74.6 

— 

'O' 

21 

-67.3 

— 

“The  electrostratic  potential  minima,  all  of  which  are  asssociated  with  the  ring  nitrogens,  were 
calculated  in  this  work.  The  p Ka  values  are  experimentally  determined,  and  were  taken  from  [  1  ] . 


that  the  ring  carbons  are  considerably  less  reactive  toward  electrophiles  than 
in  benzene,  but  more  susceptible  to  nucleophilic  attack  [1];  both  tendencies 
increase  with  the  number  of  ring  nitrogens. 

The  strongly  negative  electrostatic  potentials  produced  by  the  ring  nitrogens 
make  these  attractive  sites  for  electrophiles,  and  indeed  a  variety  of  such  re¬ 
actions  have  been  found  to  take  place  [1  ].  The  magnitudes  of  these  potentials 
decrease  with  the  number  of  nitrogens  in  the  ring,  as  shown  in  Table  2,  reflect¬ 
ing  the  fact  that  as  more  nitrogens  compete  for  the  polarizable  electronic  charge, 
each  receives  a  smaller  share  of  it.  This  phenomenon  has  previously  been  noted 
in  other  systems  containing  electron-attracting  atoms  or  groups  [12,14-16]. 
Correspondingly,  the  reactivities  toward  electrophiles  also  decrease  as  the 
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Fig.  8.  Electrostatic  potential  (in  kcal  mol-1)  in  the  plane  1.75  A  above  the  molecular  plane  of 
nitrobenzene,  XVIII.  The  projections  of  the  nuclear  positions  are  indicated.  Dashed  contours 
correspond  to  negative  values.  The  locations  of  the  most  negative  potentials  in  this  plane  are 
shown,  and  the  magnitudes  are  given  below  the  Figure.  The  planar  maximum,  above  the  C-N02 
bond  region,  is  +  14.4  kcal  mol-1. 

o 

Fig.  9.  Electrostatic  potential  (in  kcal  mol- 1 )  in  the  plane  1.75  A  above  the  molecular  plane  of  3- 
nitropyridine,  IX.  The  projections  of  the  nuclear  positions  are  indicated.  Dashed  contours  corre¬ 
spond  to  negative  values.  The  locations  of  the  most  negative  potentials  in  this  plane  are  shown, 
and  the  magnitudes  are  given  below  the  figure.  The  planar  maximum,  above  the  C-N02  bond 
region,  is  +20.9  kcal  mol-1. 
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Fig.  10.  Electrostatic  potential  (in  kcal  mol-1)  in  the  plane  1.75  A  above  the  molecular  plane  of 
4jnitropyridazine,  XI.  The  projections  of  the  nuclear  positions  are  indicated.  Dashed  contours 
correspond  to  negative  values.  The  locations  of  the  most  negative  potentials  in  this  plane  are 
shown,  and  the  magnitudes  are  given  below  the  figure.  The  planar  maximum,  above  the  C-N02 
bond  region,  is  +29.8  kcal  mol-1. 

Fig.  11.  Electrostatic  potential  (in  kcal  mol-1)  in  the  plane  1.75  A  above  the  molecular  plane  of 
2-nitropyrazine,  XII.  The  projections  of  the  nuclear  positions  are  indicated.  Dashed  contours 
correspond  to  negative  values.  The  locations  of  the  most  negative  potentials  in  this  plane  are 
shown,  and  the  magnitudes  are  given  below  the  Figure. 
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Fig.  12.  Electrostatic  potential  (in  kcal  mol-1)  in  the  plane  1.75  A  above  the  molecular  plane  of 
nitro-s-triazine,  XVI.  The  projections  of  the  nuclear  positions  are  indicated.  Dashed  contours 
correspond  to  negative  values.  The  locations  of  the  most  negative  potentials  in  this  plane  are 
shown,  and  the  magnitudes  are  given  below  the  Figure. 

number  of  ring  nitrogens  increases  [1];  for  example,  the  basicities  follow  ex¬ 
actly  the  same  trend  as  the  potential  minima,  as  seen  from  the  p Ka  values  in 
Table  2.  The  experimentally  observed  proton  affinities  of  pyridine  and  the 
three  diazines  also  decrease  in  this  order;  thus  the  observed  trend  in  this  prop¬ 
erty  can  be  explained  without  invoking  hydrogen  bonding  in  the  protonated 
form,  as  was  done  earlier  [111- 

Electrostatic  potentials  of  nitroazines 

The  substitution  of  a  nitro  group  on  an  azine  ring  introduces  another  strong 
competitor  for  the  polarizable  electronic  charge.  Electron  withdrawal  by  -N02 
is  primarily  through  induction  [2],  although  a  minor  degree  of  conjugation 
does  take  place  [17],  as  represented  below  for  the  nitropyridines. 
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It  is  to  be  anticipated,  therefore,  that  the  ring  nitrogens  will  be  less  negative 
than  in  the  unsubstituted  azines,  and  similarly  that  the  -N02  group  will  be 
less  negative  than  in  nitrobenzene,  XVIII.  Both  of  these  expectations  are  borne 
out  by  our  calculated  electrostatic  potentials. 

The  negative  lone-pair  potentials  of  the  ring  nitrogens  are  significantly 
weakened  in  going  from  the  azines  to  their  mononitro  derivatives.  This  can  be 
seen,  for  instance,  by  comparing  Figs.  5  and  6;  the  values  of  the  minima  change 
from  84.1  kcal  mol-1  in  pyrimidine,  IV,  to  —70.2  kcal  mol-1  in  5-nitropyr- 
imidine,  XV.  The  latter  value  is  fairly  typical  for  the  nitrodiazines,  except 
when  there  is  an  overlapping  and  reinforcement  of  adjacent  negative  regions, 
as  in  X  (Fig.  3 ) ,  XI  and  XIII.  For  the  nitropyridines,  VII-IX,  the  minima  are 
stronger,  approximately  —79  kcal  mol-1  (compared  to  —92.9  in  pyridine), 
while  for  XVI  and  XVII  they  range  from  —53.5  to  —63.8  kcal  mol-1.  Thus 
nitration  should  be  accompanied  by  a  general  decrease  in  the  basicities  of  the 
azines,  as  has  indeed  been  observed  for  the  nitropyridines  [1  ]. 

The  nitro  groups  in  VII-XVIII  have  their  customary  extended  regions  of 
negative  electrostatic  potential  near  the  oxygens  [12,13,16,17].  These  have 
four  minima  (Figs.  6  and  7 ) ,  which  can  be  interpreted  as  reflecting  the  oxygen 
lone  pairs.  The  magnitudes  of  these  minima  in  nitrobenzene,  XVIII,  are  —41, 
—  60,  —  60  and  —  41  kcal  mol -  however,  the  presence  of  ring  nitrogens  weak¬ 
ens  these  to  —35,  —51,  —51  and  —35  kcal  mol'7  in  4-nitropyridine,  VIII,  for 
example,  and  even  further  to  —30,  —48,  —48  and  —30  kcal  mol-1  in  5-nitro- 
pyrimidine,  XV  (Fig.  6).  Exceptions  to  these  trends  do  of  course  appear  when 
there  is  overlap  of  negative  regions,  as  in  X  (Fig.  7). 

The  electrostatic  potentials  of  all  the  nitro  derivatives,  VII-XVIII,  are  pos¬ 
itive  everywhere  above  and  below  the  ring.  This  is  true  of  nitroaromatics  even 
when  there  are  no  nitrogens  in  the  ring  [13,17].  These  positive  regions  are 
stronger  for  the  nitroazines,  however,  as  can  be  seen  by  comparing  Figs.  8  and 
9,  and  the  difference  becomes  even  more  marked  as  the  number  of  ring  nitro¬ 
gens  increases  (Figs.  10  and  11).  Thus,  while  nitroazines  should  not  undergo 
ring  interactions  with  electrophiles,  susceptibility  toward  nucleophilic  attack 
is  certainly  to  be  anticipated. 

In  all  the  nitroazines  these  positive  potentials  tend  to  concentrate  above  the 
central  regions  of  the  rings,  increasingly  so  for  nitro-s-triazine  and  nitro-s- 
tetrazine  (Figs.  9-12).  However,  in  the  nitropyridines,  VII-IX,  and  in  two  of 
the  nitrodiazines,  XI  and  XV,  the  strongest  build-up  of  positive  potential  in 
this  plane  occurs  above  the  C-N02  bond  region  (see  Figs.  9-11 ).  This  feature 
was  previously  observed  in  the  nitro  derivatives  of  a  variety  of  unsaturated 
systems,  including  aromatics  [13,17-19],  dioxins  [20],  dioxazines  and  diox- 
adiazines  [12],  acetylene  [16]  and  ethylene  [21].  (This  can  be  seen  in  Fig.  8 
for  the  case  of  nitrobenzene. )  We  have  shown  that  these  positive  build-ups  can 
act  as  initial  sites  for  nucleophilic  attack,  as  in  hydrolysis  [18].  The  magni¬ 
tudes  of  build-up  depends  very  much  upon  the  remainder  of  the  system.  Thus 
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the  electron-attracting  effect  of  one  ring  nitrogen  in  IX  causes  this  planar 
maximum  to  increase  from  14  kcal  mol- 1  in  nitrobenzene,  XVIII,  (Fig.  8)  to 
21  kcal  mol  - 1  (Fig.  9 ),  while  the  second  ring  nitrogen  in  XI  produces  a  further 
increase  to  30  kcal  mol-1  (Fig.  10).  Even  more  positive  build-up  was  seen  in 
dinitro  derivatives,  reaching  35  kcal  mol-1  in  dinitroacetylene  [16]  and  ap¬ 
proximately  50  kcal  mol-1  in  the  dinitro- 1,4-dioxadiazines  [12]. 

SUMMARY 

Our  primary  focus  in  this  paper  has  been  upon  the  relative  stabilities  and 
the  electrostatic  potentials  of  mononitroazines.  We  found  that  the  preferred 
site  for  -N02  substitution  is  on  a  carbon  /?  to  a  ring  nitrogen,  while  the  a 
position  is  the  least  stable.  For  both  the  unsubstituted  diazines  and  their  nitro 
derivatives,  however,  the  most  important  destabilizing  factor  is  the  presence 
of  two  adjacent  ring  nitrogens. 

The  introduction  of  the  strongly  electron -withdrawing  nitro  group  weakens 
the  negative  electrostatic  potentials  associated  with  the  ring  nitrogens,  which 
should  accordingly  be  significantly  less  basic  than  in  the  corresponding  unsub¬ 
stituted  azines.  (In  the  latter,  the  observed  basicities  follow  the  same  trend  as 
the  magnitudes  of  the  negative  potentials  near  the  ring  nitrogens.)  None  of 
the  nitroazines  has  any  negative  regions  above  or  below  the  ring;  instead  there 
are  positive  build-ups  above  the  central  portions  of  the  rings  and  sometimes 
above  the  C-N02  bonds,  which  are  expected  to  provide  channels  for  nucleo¬ 
philic  attack. 
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